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INTRODUCTION 


If  a  metal  contains  very  few  impurities  and  is  mechanically 
tested  in  an  inert  environment,  its  fracture  mode  will  be  either 
transgranular  cleavage  or  transgranular  ductile  microvoid  coales¬ 
cence.  If,  however,  the  metal  is  of  commercial  purity,  the  brittle 
fracture  mode  can  switch  from  cleavage  to  intergranular  separation. 
This  transition  is  undesirable  because  intergranular  fracture  is 
generally  a  very  low  energy  process. 

The  occurrence  of  intergranular  fracture  in  a  specific  situa¬ 
tion  will  depend  on  a  number  of  metallurgical  variables:  grain 
size,  yield  strength,  test  temperature,  and  notch  radius.  However, 
in  all  reported  cases,  one  observation  is  common.  Impurities  which 
have  a  very  low  solubility  in  the  bulk  have  segregated  to  the  grain 
boundaries  and  weakened  them.  The  bulk  concentration  of  these 
impurities  need  only  be  several  liundred  ppm  for  grain  boundary 
concentrations  of  5-10  atomic  percent  to  be  observed.  Examples  of 
this  impurity-induced  weakening  of  grain  boundaries  include  sulfur 
in  iron  and  nickel,  phosphorus,  tin,  and  antimony  in  steel,  bismuth 
in  copper,  and  oxygen  in  refractory  metals. 

If  chf.  cuwircnr.ient  is  not  inert,  intergranular  separation  may 
occur  even  tliough  the  fracture  mode  in  an  inert  environment  would 


be  cransgrar.uiar .  One  .-such  studied  exacitJie  of  this  is  hydrogen 
enbrictietienc .  If  s  sample  is  tasted  either  in  hydrogen  gas,  or  is 
cathodicaily  oharged  with  hydrogen  ind  then  tested  in  an  inert 
environment,  one  finds  that  the  mechanical  properties  of  the  mace- 
tial  are  often  degraded  and  that  the  fracture  mode  has  become  more 
intergranular.  As  with  impurity-induced  embrittlement,  many  addi¬ 
tional  variables  can  affect  the  aegree  to  which  hydrogen  will  cause 
embrittlement,  and  one  of  Che  worst  possible  cases  is  to  have  both 
hydrogen  and  impurities  segregated  at  Che  grain  boundaries.  However, 
it  is  clear  that  hydrogen  must  segregate  at  interfaces  and  reduce 
the  cohesive  strength  across  them 

We  believe  that  an  important  contribution  to  understanding 
some  aspects  of  these  complex  problems  of  embrittlement  in  metals 
can  be  made  through  theoretical  studies  of  bonding  in  model  systems. 
The  philosophy  is  to  isolate  groups  of  atoms  in  clusters  which 
simulate  a  possible  local  environment  of  interest  and  to  perform 
detailed  quantum  mechanical  calculations  using  molecular  orbital 
theory  to  probe  Che  nature  of  the  chemical  bonding. 

After  discussing  molecular  orbital  theory  and  the  Xa-scattered- 
wave  technique  in  the  next  section,  we  consider  in  the  following 
sections  two  recent  applications  of  this  approach  to  models  of 
potential  metallurgical  interest. 

^K)LECULAR  ORBITAL  THEORY  AND  THE  SCF-Xa-SW  CLUSTER  METHOD 

The  most  characteristic  feature  of  molecular  orbital  theory  is 
that  it  provides  an  independent-particle  interpretation  of  the 
electronic  structure  of  a  molecule  or  solid.  That  is,  each  electron 
of  the  system  has  associated  with  it  a  characteristic  spin-orbital 
and  energy  level  which  are  determined  by  solving  Schrodinger 's 
equation  for  an  electron  in  the  average  potential  energy  of  all  the 
other  electrons  in  the  system.  The  most  familiar  of  such  techniques 
is  the  Hartree-Fock  self-consistent  field  method,  in  which  each 
molecular  orbital  is  traditionally  represented  by  a  linear  combina¬ 
tion  of  atomic  orbitals  (the  HF-SCF-LCAO  method).  For  all  but  the 
simplest  molecules  one  must,  at  least  for  the  present,  choose 
atomic-orbital  basis  sets  which  are  fairly  far  from  the  correct 
solutions  of  the  Hartree-Fock  Schrodinger  equations.  We  refer  to 
this  type  of  approximate  solution  of  the  Hartree-Fock  equations 
using  the  LCAO  representation  as  the  AHF-SCF-LCAO  method. 

It  is  well  known  that  the  Hartree-Fock  theory  provides  a  poor 
description  of  the  electronic  structure  of  a  molecule  or  solid  at 
large  internuclear  distances  and  that  it  leads  to  an  improper 
description  of  tiie  dissociation  of  such  systems.  In  order  to  cor¬ 
rect  for  this  behavior,  it  is  necessary  for  one  to  work  with  linear 
combinations  of  Slater  determinants  of  molecular  orbitals,  the 


so-called  "conf iyuracion-incaraccion"^  (CI'i  procedure,  vhich  is  cor 
only  conipiicacea  compucactonally ,  out  also  iescroys  rhe  indepencccr 
electron  incerprecaricn .  A  furtner  diificulry  is  that  in  perrorruLr. 
numerical  calculacions  using  che  LCaO  approximacion  at  either  the 
AHF  or  Cl  levels  of  sophistication,  the  nunber  of  rvo-electron 
multicenter  integrals  which  must  be  calculated  even  for  moderately 
sized  systems  becomes  enormous.  For  example,  in  a  recent  calcula¬ 
tion  on  a  molecule  consisting  of  26  atoms,  ail  of  atomic  numbers 
less  Chan  ten,  10^-  integrals  had  to  be  evaluated  and  the  calcula¬ 
tions  required  the  equivalent  of  192  h  on  an  IBM  360/195  computer , 
one  of  the  fastest  machines  available. 

A  much  more  practical  independent-electron  theory,  which  avoid 
most  of  the  difficulties  associated  with  the  AHF-SCF-LCAO  and  Cl 
methods,  is  che  recently  developed  SCF-I\a-SW  cluster  method.  Based 
on  the  combined  use  of  Slater's  Xa  statistical  theory  of  exchange 
correlation  and  Johnson's  multiple-scattered-wave  formalism,  this 
technique  allows  one  to  calculate  accurate  self-consistent  spin- 
orbital  energies,  wave  functions,  charge  densities,  and  total 
energies  for  comple.x  molecules  and  clusters,  yet  requires  only 
relatively  small  amounts  of  computer  time.  Like  the  AHF-SCF-LCAO 
approach,  che  SCF-Xa-SW  method  is  an  approximate  theory  of  many- 
electron  systems,  but  it  is  not  an  approximation  to  Hartree-Fock 
theory.  First  of  all,  in  contrast  to  KF  theory,  the  Xa  statistical 
expression  for  the  total  many-electron  energy  automatically  goes 
to  the  proper  separaced-acom  limit  as  the  internuclear  distances 
are  increased  to  infinity,  i.e.,  as  the  system  dissociates.  Second 
it  has  been  shown  chat  the  Xa  theory  satisfies  Fermi  statistics, 
thereby  ensuring  the  proper  ordering  and  occupation  of  electronic 
energy  levels,  while  the  HF  theory  does  not.  Third,  it  has  been 
proven  that,  under  che  proper  conditions,  the  Xa  statistical  total 
energy  rigorously  satisfies  both  the  virial  and  Hellmann-Feynman 
theorems,  which  facilitates  the  calculation  of  the  equilibrium 
cohesive  energies  of  molecules  and  solids.  Finally,  the  SCF-Xa-SV 
method,  in  conjunction  with  the  "transition-state"  concept,  leads 
to  an  accurate  approximate  description  of  the  excited  electronic 
states  and  charge  distributions  of  polyatomic  systems,  including 
the  effects  of  orbital  relaxation  which  are  usually  neglected  in 
applying  Koopmans  theorem  in  the  AHF-SCF-LCAO  approach. 

The  SCF-Xa-SW  method  and  transition-state  procedure  have 
already  been  applied  successfully  to  a  wide  range  of  molecules  and 
clusters  in  solids.  Because  the  required  computational  effort  does 
not  increase  Inordinately  with  the  number  of  electrons  per  atom  or 
with  the  number  of  atoms  in  the  molecule  or  cluster,  the  SCF-Xi-SW 
method  is  ideally  suited  for  describing  the  electronic  structures 
of  transition-metal  complexes  and  clusters  representing  the  local 
mo]ecular  environment  of  a  grain  boundary. 


MCLECULA?.  OREir.iL  ^^ODEL  FOR  GR.\.IN 
3Y  SULFUR 


EM3Rn 


■JF  NICKEL 


As  discussed  in  Che  last  section,  the  theoretical  lethod  we 
use  to  investigate  the  electronic  structure  of  the  local  environ¬ 
ment  at  a  grain  boundary  consists  of  two  parts.  The  first  is  to 
use  a  cluster  of  atoms  to  represent  the  local  environment  and  the 
second  is  to  use  molecular  orbital  theory  to  solve  for  the  elec¬ 
tronic  structure  of  the  cluster  of  atoms. 


To  apply  this  technique  to  grain  boundaries  one  must  first 
choose  a  cluster  that  resembles  in  some  way  Che  local  environment 
an  impurity  atom  might  have  were  it  at  the  grain  boundary.  The 
choice  of  such  clusters  is  considerably  simplified  due  to  the 
recent  work  of  Ashby,  Spaepen,  and  Williams,  who  have  shown  chat 
the  structure  of  grain  boundaries  in  fee  solids  can  be  described  in 
terras  of  atoms  at  Che  vertices  of  Bernal  deltahedra.  The  individual 
deltahedra  are  small  clusters  containing  four  to  thirteen  atoms  and 
thus  provide  a  natural  initial  model  of  local  atomic  structure  with 
which  to  investigate  the  electronic  properties  by  SCF-.Ya-SW  calcu¬ 
lations.  In  order  to  investigate  the  feasibility  of  applying  such 
methods  to  this  problem  we  have  chosen  to  use  one  particular  cluster 
in  this  work. 


The  particular  cluster  we  chose  is  the  tetragonal  dodecahedron 
shown  in  Fig.  1.  It  contains  eight  atoms  and  has  Dgd  point  group 
symmetry.  We  do  not  claim  that  this  structure  is  necessarily  a 
prevalent  one  in  grain  boundaries.  However,  it  does  contain  four 
atoms  (numbered  1,  2,  3,  and  4)  slightly  but  symmetrically  displaced 
from  the  x-y  plane  and  fd’ur  other  atoms  further  displaced  from  this 
plane.  One  can  roughly  consider  that  the  four  atoms  nearest  the 
x-y  plane  are  grain  boundary  atoms  and  that  Che  other  four  are  in 
the  first  atomic  layers  away  from  the  precise  boundary. 

It  is  well  established  experimentally  that  sulfur  embrittles 
nickel.  In  order  to  investigate  the  electronic  effects  of  a  sulfur 
atom  at  a  grain  boundary  in  Ni,  we  have  carried  out  two  calculations 
on  a  model  system  in  which  the  Ni  atoms  were  placed  at  each  of  the 
eight  sites  of  the  cluster  shown  in  Fig.  1  with  an  internuclear 
separation  of  2*49  A.  In  the  first  calculation,  we  considered  the 
cluster  containing  only  the  eight  nickel  atoms.  For  the  second 
calculation,  a  sulfur  atom  was  placed  at  the  center  of  the  cluster. 
An  interstitial  site  of  this  type  is  a  possible  position  which  a 
sulfur  atom  would  choose  at  a  grain  boundary.  For  these  two  con¬ 
figurations  we  have  determined  the  molecular  orbitals  for  the 
clusters,  their  corresponding  energy  levels,  and  the  valence  charge 
densities . 

The  energy  levels  for  i^ach  cluster  arc  shown  in  Fig.  2.  We 
note  that  in  the  NigS  cluster  several  of  the  energy  levels  have 
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been  lowered  subs rancially .  Those  chac  have  been  lowered  have 
orbital  wave  functions  which  are  highly  localized  on  the  sulfur 
atoci  and  its  four  nearest  neighbors,  i.e.,  the  four  nickel  stocos 
closest  to  the  x-y  plane.  3ased  on  this  sixpie  analysis  one  night 
assuEe  that  soxe  Mi-Hi  bonds  have  been  weakened. 

The  bonding  nay  be  better  understood  if  one  analyzes  the 
molecular  orbitals  of  the  clusters  in  more  detail;  however,  in 
this  particular  case  a  sore  succinct  picture  of  the  bonding  in  the 
cluster  can  be  obtained  by  sununing  Che  orbital  charge  densities 
over  all  the  occupied  valence  orbitals.  One  can  then  plot  contours 
of  the  total  valence  charge  density  which  exists  on  any  plane. 

These  are  shown  in  Fig.  3  for  the  x-z  plane.  In  this  figure  we 
have  plotted  Che  same  sec  of  contours  for  both  clusters.  The  three 
outermost  contours  (numbered  1,  2,  and  3)  in  both  clusters  surround 
all  atoms.  However,  contour  number  4  (which  is  shown  with  a 
heavier  line  for  clarity)  is  quite  different  in  the  two  clusters. 

In  Nig)  it  shows  bonding  between  the  pairs  of  nickel  atoms  1  and  7, 

7  and  8,  and  4  and  8,  where  the  characteristic  "neck"  regions 
betw’een  the  pairs  of  atoms  is  responsible  for  bonding.  In  NigS, 
the  shape  of  contour  number  4  is  quite  different.  The  character¬ 
istic  neck  bonding  region  between  atoms  1  and  7,  and  4  and  8  is 
gone  and  consequently  so  is  the  concomitant  strong  bonding  between 
the  two  types  of  nickel  atoms.  Furthermore,  the  size  of  the  neck 
region  between  atoms  7  and  8  is  reduced  from  that  in  Nig,  implying 
a  weakening  of  this  metal-metal  bond  as  well.  Contour  4  in  NigS 
is  split  into  two  distinct  regions.  The  upper  region  of  contour  4 
and  an  inner  contour  of  the  nickel  atoms  nearest  the  sulfur  atom 
enclose  both  the  nickel  atoms  and  the  sulfur.  This  would  lead  to 
very  strong  bonds  between  this  set  of  nickel  atoms  and  the  sulfur. 
However,  there  is  no  such  bonding  interaction  berween  the  sulfur 
atoms  and  the  other  nickel  atoms. 

From  these  results  a  clear  picture  begins  to  emerge.  Adding 
sulfur  to  the  cluster  of  nickel  atoms  causes  strong  NiS  bonds  to 
be  formed  between  the  nickel  atoms  nearest  the  sulfur.  Yet  this 
weakens  bonds  which  the  nickel  atoms  would  ordinarily  have  with 
the  other  nickel  atoms  that  are  not  as  close  to  the  sulfur  atom. 

In  this  cluster  the  strong  Ni-S  bonds  involve  the  sulfur  atom  and 
the  four  nickel  atoms  nearest  the  x-y  plane.  We  suggested  earlier 
that  these  atoms  might  be  thought  of  as  being  in  the  grain  boundary. 
However,  the  bonds  between  the  nickel  atoms  nearest  the  .x-y  plane 
and  those  farther  away  are  weakened.  Those  bonds  might  represent 
those  from  the  boundary  out  to  the  first  layer  of  the  bulk.  Clearly, 
if  those  are  weakened,  the  stress  required  for  fracture  would  be 
less . 


Therefore,  based  on  the  calculations  wc  can  suggest  an  elec¬ 
tronic  mechanism  by  which  grain  boundary  ombri  ttlcnicnt  could  occur. 
We  stress  that  these  are  the  results  from  calculations  on  only  one 


sysceni.  rurchertroria ,  ve  nave  chosen  an  embrlrrling  aien:enc  vhicn 
probably  occupies  an  incersciPiai  size  ar  one  grain  boundair.’. 

Other  eruD ric tiin 2  elenenrs  vould  ciaariy  be  subs ticucionai .  Also, 
many  other  types  of  geocetricai  conf iguracicn  are  possible  at  tne 
boundary/.  Vet,  with  these  results,  we  believe  that  it  is  not  pre¬ 
sumptuous  to  propose  that  the  above-described  mechanism  is  :ne 
mechanism  by  which  segregated  elements  embrittle  the  grain  bounder;/. 

.VOLECUL.AR  ORBITAL  MODEL  FOR  HYDROGE"  EMERITTLEIE.N’T  OF  IRON 

The  causes  and  mechanisms  of  hydrogen  embrittlement  have  been 
and  will  continue  to  be  the  subjects  of  detailed  investigation. 

As  a  result,  several  models  have  been  offered  in  explanation  of  the 
embrittlement  phenomenon.  These  models,  however,  cannot  account 
for  embrittlement  at  the  atomic  level,  as  a  chemical  interaction 
between  hydrogen  and  the  embrittled  metal.  It  is  the  purpose  of 
this  paper  to  present  the  results  of  SCF-Xa-SW  molecular-orbital 
calculations  on  iron-hydrogen  clusters  in  order  to  e.xplain  hydrogen 
embrittlement  of  ferrous  alloys  at  the  most  fundamental  level. 

Hydrogen  atoms  are  most  frequently  found  in  octahedral  holes 
of  transition  metals.  If  molecular  hydrogen  were  adsorbed  within 
tetrahedral  holes  then  concerted  dissociation  of  molecular  hydrogen 
would  result  in  placing  atomic  hydrogen  in  octahedral  holes. 
Accordingly,  the  four-atom  iron  cluster  shown  in  Fig.  4  will  be  -used 
to  model  the  tetrahedral  hole  and  the  site  for  hydrogen  adsorption 
in  bcc  iron.  (Both  the  tetrahedral  and  octahedral  holes  of  bcc 
iron  are  misnamed  since  they  actually  present  a  D2d  and  tetragonal 
symmetry  respectively.)  In  an  earlier  paper,  the  results  of  SCF- 
Xa-SW  molecular-orbital  calculations  of  four-,  nine-,  and  fifteen- 
atom  clusters  of  iron  were  presented.  It  was  shown  that  these 
clusters  provide  a  model  for  the  local  electronic  structure  of  bcc 
iron.  .A  comparison  of  these  results  indicates  that  both  the  four- 
atom  and  fifteen-atom  clusters  exhibit  the  large  exchange  splitting 
and  high  magnetic  moment  characteristic  of  bulk  iron.  Net  spin 
electron  number  per  atom  for  Fe]_5  is  2.7  and  for  the  Fe^  cluster 
2.5  compared  with  2.12  for  the  bulk.  The  density  of  states  in  the 
four-atom  cluster  is  sufficiently  similar  to  the  density  of  states 
in  the  15-atom  cluster  to  justify  the  use  of  Fe4  as  a  model  for  a 
site  of  hydrogen  adsorption.  It  is  important  to  note  that  the 
fifteen-atom  cluster  provides  three  environments  for  the  atoms  of 
the  cluster.  In  bulk  iron  (neglecting  surface  effects)  each  iron 
atom  is  in  the  same  environment.  It  is  this  fact  which  is  respon¬ 
sible  for  the  variation  of  properties  predicted  by  the  cluster 
method  with  those  measured  in  bulk  iron.  In  a  calculation  of  this 
type,  however,  where  an  impurity  is  being  modeled,  the  im.purity 
breaks  the  symmetry  of  the  bulk  and  there  are  a  number  of  environ¬ 
ments  around  the  impurity  site.  As  a  result,  the  cluster  method  is 
ideally  suited  for  the  study  of  dilute  impurity  states  i,id 
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The  nolecular-orbital  energy  levels  or  a  hydrogen  xolecuie 
inserted  into  the  four-acora  cluster  is  shoun  in  Fig.  3.  The  aj_ 
and  h?  oroitals  of  the  four-acor:!  iron  cluster  have  th.e  sat'.e  svti- 
metry  as  the  hydrogen  nolecuie  bonding  and  antibonding  orbitals. 
These  hydrogen  molecular  orbitals  couple  with  the  lowest  lying 
ai  and  b2  orbitals  of  the  Fe4  cluster  to  produce  the  molecular 
orbitals  shown  in  Figs.  6  and  7.  As  a  result  of  this  coupling  the 
lowest  lying  and  bn  orbitals  of  Fa^Hi  are  stabilized  relative 
to  the  position  of  these  orbitals  in  Fe^.  The  a^  orbital  is  pre¬ 
dominantly  a  F.-H  bonding  orbital  the  H-Fe  interaction  is  non- 
bonding.  The  b2  orbital  in  the  Fe^  cluster  is  d-d  --bonding  between 
nearest  neighbors  and  o-bonding  between  second-nearest  neighbors. 

The  addition  of  H2  strengthens  the  bonding  interaction  between 
second-nearest  neighbors,  has  little  effect  on  the  nearest-neighbor 
interaction,  and  is  strongly  antibonding  between  the  hydrogen  atoms. 
The  H-H  distance  in  this  cluster  is  the  molecular  H-H  distance. 

This  is  the  distance  which  is  realized  when  there  is  no  electron 
density  in  a  H-H  antibonding  orbital.  In  this  cluster  some  of  the 
electron  density  from  the  metal  cluster  b2  orbital  is  transferred 
into  the  antibonding  orbital  of  the  hydrogen  molecule.  Tlie  hydrogen 
molecule  H-H  distance  is  no  longer  a  stable  equilibrium  distance 
and  the  hydrogen  atoms  will  begin  to  dissociate.  All  transition 
metals  are  known  to  dissociate  hydrogen  molecules.  The  efficiency 
of  this  dissociation  will  depend  on  the  amount  of  overlap  between 
the  metal  orbital  of  suitable  symmetry  and  the  antibonding  orbital 
of  the  hydrogen  molecule.  Other  than  the  low-lying  a^^  and  b2 
orbitals,  no  other  orbitals  of  this  cluster  have  any  significant 
hydrogen  admi.tture.  Two  additional  spin  orbitals  are  occupied, 
however,  by  the  two  electrons  contributed  by  the  hydrogen  molecule. 

Figure  8  gives  the  molecular  orbital  energy  levels  for  the 
same  arrangement  of  iron  atoms  but  now  the  H-H  distance  has  been 
increased  as  a  consequence  of  the  electron  density  in  the  H-H 
antibonding  orbital.  Again  the  lowest  lying  a^  orbital  is  the  only 
orbital  with  significant  H-H  bonding  contributions  (Fig,  9).  The 
lowest  lying  b2  orbital  is  H-H  antibonding,  as  in  the  previous 
case,  but  this  interaction  has  been  lessened  as  is  clearly  seen  in 
Fig.  10.  There  is  a  reasonably  strong  H-Fe  bonding  interaction 
within  this  b2  orbital,  and  this  interaction  should  be  maximized 
when  the  hydrogen  atom  is  sitting  between  the  iron  second-nearest 
neighbors,  i.e,,  when  the  hydrogen  atoms  are  in  the  tetragonal 
holes  of  the  bcc  structure. 

While  hydrogen  adsorption  and  dissociation  may  be  explained 
in  a  fairly  straightfoi-ward  manner  by  molecular-orbital  theory, 
attempts  Co  explain  embrittlement  at  Che  molecular  level  must  be 
assisted  by  tu'isidering  embrittlemcnr  as  a  competi  tive  process 
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between  two  nechanisns  .  Tlie  first  of  tr.ese  is  the  process  or 
plastic  flow  facilitated  by  iis location  movenenc .  The  second  is 
charactericed  by  nuclear  notion  in  a  airection  which  will  counter 
the  strain  enercy  of  the  crystal,  a  purely  elastic  resocnse.  An 
energy  of  activation  itay  be  associated  with  each  of  these  processes 
The  activation  energy  for  the  first  process  will  deper.c  preconi- 
nately  upon  the  structure  of  the  crystal,  the  Burgers  vector,  etc., 
while  the  activation  energy  of  the  second  process  is  sinpiy  the 
fracture  energy  of  the  crystal.  A  ductile  r.etal  responds  to 
excessive  strain  through  the  first  nechanism,  and  this  nay  be 
attributed  to  the  fact  that  the  activation  energy  for  plastic  flow 
is  significantly  less  than  the  activation  energy  for  brittle  frac¬ 
ture.  A  brittle  natarial  will  respond  to  strain  by  fracturing 
before  plastic  flow  begins  because  the  activation  energy  for  dis¬ 
location  movement  is  higher  than  the  activation  energy  for  brittle 
fracture.  A  ductile-brittle  transition  must  be  the  result  of  a 
mechanism  that  either  raises  the  activation  energy  for  plastic  flow 
or  lowers  the  energy  for  brittle  fracture. 

If  the  process  of  hydrogen  adsorption  affects  the  activation 
energy  for  dislocation  flow,  then  hydrogen  adsorption  must  alter 
the  structure  of  bcc  iron.  Structural  changes  are  often  manifest 
as  changes  in  the  molecular  orbital  energy  distribution  near  the 
Fermi  level.  Figure  11  shows  the  molecular  orbital  energy  distribu 
tion  near  the  Fermi  level  for  a  tetrahedron  of  iron  atoms,  the 
four-atom  cluster  modeling  bcc  iron,  and  three-hydrogenated-iron 
cluster  corresponding  to  adsorbed  molecular  hydrogen,  adsorbed 
dissociated  molecular  hydrogen  and  a  hydrogen  saturated  cluster, 
Fe4H(3,  In  the  tetrahedral  clu.ster,  there  is  an  unoccupied  t2 
orbital  near  the  Fermi  level,  this  orbital  is  split  when  the  sym¬ 
metry  of  the  cluster  is  reduced  from  tetrahedral  to  Dgd  as  in  the 
A-atom  bcc  cluster.  The  e  orbital  is  stabilized  and  becomes  occu¬ 
pied  while  the  b2  orbital  is  destabilized.  This  process  is  remi¬ 
niscent  of  the  Jahn-Teller  effect  of  inorganic  complexes,  where  a 
complex  spontaneously  distorts  in  order  to  lower  its  total  energy. 
This  leads  to  the  conjecture  that  the  bcc  structure  of  iron  is  a 
consequence  of  a  Jahn-Teller  distortion  of  a  fee  structure.  When 
hydrogen  molecules  are  added  to  the  bcc  cluster  the  e  orbital  is 
depopulated  and  is  nearly  degenerate  with  the  bg  orbital,  and  is 
similar  to  the  energy  distribution  of  the  tetra'.iedral  iron  cluster. 
The  dissociation  of  the  hydrogen  molecule  stabilizes  the  b2  orbital 
but  the  addition  of  more  hydrogen  seems  once  again  to  depopulate 
this  orbital.  Hydrogen  adsorption  would  appear  to  promote  a  close- 
packed  phase  transition  near  the  site  of  adsorption.  There  is 
evidence  that  hydrogen  adsorption  does  promote  ptiase  transforma¬ 
tions.  In  HCP  iron  hiydrogen  adsorption  is  known  to  lower  tiie 
stacking  fa,ult  energy  indicating  chat  hydrogen  is  promoting  the 
fcc-hcp  phase  transition. 

This  is  just  the  t^-pe  of  effect  th  it  would  r.nise  the 
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accivaci^-'n  ener5y  for  disLocacion  move  nerve  around  the  sires  of 
fiydrogerv  adsorpeion.  However,  no  surface  recons erucr ion  of  irsn 
due  ro  hydrogen  adsorpeion  has  ever  been  noced  and  fee  merals, 
nccably  nickel,  embrictie  in  a  hydrogen  environment;,  Vnile  hvdregen 
undoubcedly  has  an  effecc  on  che  structure  of  bcc  iron,  it  is  unlike 
ly  that  this  alone  can  account  for  hydrogen  embrittlement. 

This  leaves  the  possibility  that  hydrogen  adsorpeion  lowers 
che  activation  energy  of  brittle  fracture.  If  this  is  che  case 
Chen  hydrogen-adsorbed  clusters  should  have  molecular  orbitals  that 
show  similarities  to  metal  clusters  that  have  been  strained,  for 
such  clusters  require  less  additional  energy  to  fracture,  i.e.,  che 
activation  energy  for  brittle  fracture  has  been  lowered,  .^ccordingl 
Che  molecular  orbital  energy  levels  for  che  four-atom  bcc  cluster 
in  which  che  iron  atoms  have  been  strained  20%  along  the  twofold 
rotation  axis  was  calculated.  The  relative  energy  distribution  of 
the  majority  spin  states  for  this  cluster  is  shown  in  Fig.  12, 
along  v;ith  the  four-atom  model  of  bcc  iron  and  the  dissociated 
molecular  hydrogen-iron  clusters.  (The  minority  spin  states  for 
all  clusters  show  che  same  ordering  of  the  energy  levels  as  do  the 
majority  spin  states.  Mo  information  is  lost  by  studying  only  che 
majority  spin  states,  while  the  forest  of  levels  is  considerably 
thinned , ) 

The  process  of  straining  the  bulk  material  results  in  con¬ 
siderable  rearrangement  of  the  energy  levels  within  the  energy’ 
manifold.  A  rough  measure  of  Che  difference  between  two  energy 
manifolds  is  Che  number  of  intersections  between  the  connecting 
lines  of  Che  manifolds.  There  are  ten  intersections  between  che 
bulk  and  che  strained  material,  while  there  are  only  two  intersec¬ 
tions  connecting  the  strained  and  che  hydrogenated-iron  energy 
manifolds.  Saturation  of  the  Fe^  cluster  with  hydrogen  results  in 
almost  no  rearrangement  of  the  energy  manifold.  What  rearrangement 
occurs  is  confined  to  the  Lowest  five  orbitals  of  the  manifold. 

The  direction  in  which  the  molecular-orbital  energy  levels  move 
upon  straining  can  be  seen  to  be  the  result  of  removing  electron 
density  from  between  nearest  neighbors,  and  replacing  it  in  the 
region  perpendicular  to  the  direction  of  strain,  i.e.,  between 
second-nearest  neighbors.  The  result  should  be  a  contraction  of 
the  cluster  perpendicular  to  che  direction  of  strain.  This  is  the 
Poisson  effect .  Within  this  analysis  it  is  clear  that  hydrogen  atoms 
in  tetragonal  holes  will  mimic  the  behavior  of  strain  by  increasing 
the  electron  density  between  second-nearest  neighbors.  It  is  worth 
noting  chat  only  che  two  lowest  lying  orbitals  of  the  hydrogenated- 
iron  cluster  have  anv  hydrogen  admixture;  therefore,  the  rearrange¬ 
ment  of  the  other  orbitals  must  be  the  result  of  a  relaxation  effect 
due  to  the  stabilization  of  che  lowest  and  b,,  orbitals. 

It  seems  reasonable  then  that  if  the  lowest  lying  aj^  and  '02 


orbitali  of  the  scrair.ed  and  hvdrcgenarei  oLuscer  are  similar  im 
cocpositioa,  tiie  ocner  orbitals  of  the  maaifold  should  also  sm:v 
similari ties  .  ri’ures  13  and  1-  acmoara  r;'e  lowest  bi  arbitals  Ir. 

an  xz  plane  containing  the  second  neighbors  tor  the  uns trainee  :na 
strained  (20'"^^  :e^  clusters,  resneccively .  Again,  bocn  the  h'-’-iroge 
adsorbed  (Fig.  10)  and  strained  cluster  'Fig.  li)  show  greater 
electron  density  between  second  neighbors  ana  even  similar  e_ectrjn 
distributions  in  areas  away  from  the  hydrogen  atom. 
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